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Cosmic microwave background (CMB) polarization observations will require superb control of 
systematic errors in order to achieve their full scientific potential, particularly in the case of at- 
tempts to detect the B modes that may provide a window on infiation. Interferometry may be a 
promising way to achieve these goals. This paper presents a formalism for characterizing the effects 
of a variety of systematic errors on interferometric CMB polarization observations, with particular 
emphasis on estimates of the B-mode power spectrum. The most severe errors are those that couple 
the temperature anisotropy signal to polarization; such errors include cross-talk within detectors, 
misalignment of polarizers, and cross-polarization. In a _B mode experiment, the next most serious 
category of errors are those that mix E and B modes, such as gain fluctuations, pointing errors, and 
beam shape errors. The paper also indicates which sources of error may cause circular polarization 
(e.g., from foregrounds) to contaminate the cosmologically interesting linear polarization channels, 
and conversely whether monitoring of the circular polarization channels may yield useful informa- 
tion about the errors themselves. For all the sources of error considered, estimates of the level of 
control that will be required for both E and B mode experiments are provided. Both experiments 
that interfere linear polarizations and those that interfere circular polarizations are considered. The 
fact that circular experiments simultaneously measure both linear polarization Stokes parameters 
in each baseline mitigates some sources of error. 

PACS numbers: 95.75.Hi,95.75.Kk,95.85.Bh,95.85.Fm,98.70.Vc,98.80.-k,98.80.Es 



INTRODUCTION 



Cosmic microwave background (CMB) polarimetry is 
one of the most exciting frontiers in cosmology. CMB 
polarization has already been detected Q, |M 13 ^ IM IS i 
and we may expect future instruments to characterize 
the polarization signal in much greater detail (e.g., 0)- 
In the near future, CMB polarization data are expected 
to refine estimates of cosmological parameters Q , probe 
the ionization history of the Universe and the details 
of recombination \l(J^ , and measure gravitational lensing 
due to large-scale structure ^J. Most exciting of all, 
polarization maps may provide a direct probe of an infla- 
tionary epoch in the extremely early Universe by detect- 
ing the sig nature of primordial gravitational radiation 

laiiiSiis 

A crucial insight into the analysis of CMB polariza- 
tion data is the fact that any CMB polarization map 
can be divided into two components, a scalar component, 
traditionally denoted E, and a pseudoscalar component 
called B. The CMB is weakly polarized, meaning that 
both of these components are much smaller than the un- 
polarized (temperature) anisotropy. Furthermore, the B 
component is expected to be much weaker than E, since 
scalar densit y pe rturbations produce only E to linear or- 
der [Hill "El m. (See Figure □) Experiments to 
date have detected only the E component. In the future, 
the search for the weaker _B-type polarization will be a 
high priority, as the B modes may contain the imprint of 
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gravitational waves produced during inflation. 

Characterization of CMB polarization requires both 
very low noise and exquisite control of systematic errors. 
In particular, some sources of systematic error may cause 
the polarization signal to be contaminated by the much 
larger unpolarized anisotropy, while others mix the E 
and B components. As efforts to design i?-mode experi- 
ments intensify, it is important to consider carefully the 
susceptibility of different designs to various kinds of er- 
ror. Hu et al. ^3 have provided a detailed framework for 
performing such an analysis in the context of an imag- 
ing experiment. For interferometric measurements, the 
issues are somewhat different. The purpose of this paper 
is to forecast the effects of a variety of systematic errors 
on interferometric measurements. 

Interferometric methods have played an important 
role in measurements of CMB anisotropy and polariza- 
tion. Pioneering attempts to detect CMB anisotropy 
with interferometers are described in ^3 and ^3- Sev- 
eral groups have successfully detected primary CMB 
anisotropics jlO, HI il IHlil and polarization 0, |1| 
using interferometers. The formalism for analyzing CMB 
data from interferometers has been developed by a num- 
ber of authors IH IH Ulil IH Is^l as well as in the 
experimental papers cited above. 

In any data set that fails to cover the entire sky, it is 
impossible to separate the E and B components perfectly 
|3ll Isa . 133 . l34| . The operation of separating a polariza- 
tion map into E and B components is nonlocal when 
the map is viewed in real space, but in Fourier space or 
spherical harmonic space, it can be done locally (mode 
by mode). Since interferometric data sample the sky in 
the Fourier domain, E-B separation may be cleaner for 
interferometric data than for maps made with single-dish 
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FIG. 1: Power spectra for temperature anisotropy (T), TE cross-correlation (X, absolute value plotted), _E-type polarization, 
and B-type polarization. The best-fit parameters from the three-year WMAP data were used |lo] ^ 
T/S — 0.01. The right panel shows the ratios of the power spectra. 



with a tensor-to-scalar ratio 



instruments |35ll3 

As we will see, a variety of systematic errors in inter- 
ferometers can be modeled via Jones matrices 17. 37.. .3^ 
and by deviations of the antenna patterns (including 
cross-polar contributions) from assumed ideal forms. We 
will assume that each of these errors can be characterized 
by small unknown parameters, such as gain fluctuations, 
cross-talk between detectors, pointing errors, etc. We 
will first calculate the effect of each error on the measured 
visibilities. We will then provide a method of quantify- 
ing the effects of each of these errors on estimates of the 
polarization power spectra Cf, that can be obtained 
from a hypothetical data set. 

This paper has the following structure. Section ^ 
presents the mathematical formalism we will use to de- 
scribe interferometric visibilities for polarization data. 
Section IIIII presents the effects of various systematic er- 
rors on the visibilities extracted from a hypothetical 
CMB experiment. Section Hvl presents a method of fore- 
casting errors on power spectrum estimates from errors 
on visibilities. Sections IVI and IVTI contain results showing 
how the error forecasts on both E and B power spectra 
depend on the parameters that characterize the various 
systematic errors. Section IV ill presents a discussion of 
the implications of these results, and a brief appendix 
contains a useful mathematical result. 

Sections llVl through lVIl contain quite a bit of technical 
detail. The particularly busy or impatient reader should 
note that the key ideas of Section IIVI are summarized at 
the beginning, and the final results of Sections Ivl and IVTI 
are summarized in Section IVIII and Table U 



tion f . The electric field of the wave is 



E(x,i;f) =ein(f)e^(''-^-'^*), 



(2.1) 



where the wave vector k = —{uj/c)r. (Of course, as usual 
the physical electric field is just the real part of this com- 
plex quantity.) 

The complex vector ei„(f) is perpendicular to f. Its 
direction fiuctuates rapidly in time (except in the case of 
completely polarized radiation) . As usual, all observables 
will be averages taken over a time that is long compared 
to those fluctuations. To be specific, let (einXi^mF) be 
Cartesian components of ein in the plane perpendicular 
to f , and define a 2 x 2 matrix S with components 



(2.2) 



where i,j range over X,Y. The matrix S can be ex- 
pressed in terms of the standard Stokes parameters 



S = 



I + Q U + iV 
U-iV I~Q 



(2.3) 



The constant of proportionality in equation H2.2|l depends 
on the system of units being used. We will follow the 
common practice in CMB studies of expressing all Stokes 
parameters in dimensionless "AT/T" form. 

In general the electric field at any given point will be 
a superposition of waves coming in from all directions r. 
We assume that the temporal fluctuations in the incom- 
ing waves from two distinct directions are uncorrelated 
with each other: 



II. FORMALISM 



Antenna Patterns 



Consider a monochromatic plane wave of angular fre- 
quency Lo = 2'Kc/X approaching the origin from a direc- 



eini(ri)einj(f2)*) = Sij{ri)6{ii - h). 



(2.4) 



with Sij related to the Stokes parameters as above. 

Consider a single-mode antenna that is designed to 
be sensitive to only one polarization direction. We can 
model this antenna as a device that sums all of the incom- 
ing radiation, weighted by some vector-valued antenna 
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pattern a, to produce an output that is the real part of 
eout = J dHsL{r) ■ ein(f)e^(*' «--*\ (2.5) 

where ^ is the location of the antenna. 

Throughout this paper, we will consider experiments 
in which the beam width is small enough that the flat- 
sky approximation is appropriate in analyzing any sin- 
gle pointing of the instrument. (Mosaicking of multiple 
pointings of such an instrument is considered in (301].') In 
that case, we can represent the direction f by a vector 
in the plane (specifically, the tangent plane to the sphere 
at the pointing center) with components {x,y). In this 
approximation, for instance, an antenna with a Gaussian 
beam pattern that is sensitive only to linear polarization 
the X direction would have a c>c xexp[— (x^ -I- 2/^)/(4tT^)], 
while an antenna that is sensitive to either right or left 
circular polarization would have x±iy in place of x. The 
factor 4 in the Gaussian is present because we would like 
a to be the Gaussian width of |ap, not a. 

In some experiments, the same antenna may be used 
to measure two polarization states (typically either "hor- 
izontal" and "vertical" linear polarization, or left and 
right circular polarization). In that case, the output of 
the antenna will be a two-component vector Cout, and the 
antenna pattern will be a 2 x 2 matrix A: 



J (frA{r) 



i(k-|-ut) 



(2.6) 



In general, we will model antennas as sensitive to two 
polarization states in this manner. If in a particular ex- 
periment only one output is actually measured, we can 
simply ignore the other component of Cout- 

We can of course express the components of the vectors 
Cin and eout in any basis we like. In particular, we can 
resolve these vectors in either a linear polarization basis 
with components (ex, ey) or a right- and left-circular ba- 
sis with components (e_R,eL). The two bases are related 
by a unitary transformation 



with 



J_ A z 
V2 U 



(2.7) 



(2.8) 



In either case, an ideal antenna, i.e., one with equal re- 
sponse to both polarization states and no mixing between 
them, would have A equal to a scalar function times the 
identity matrix. 

In the circular polarization basis, the elements of the 
Stokes parameter matrix S become 



Scire — R-circ ' ^\in ' -E^circ (2-^) 

with Slin given by equation l|2.3|) . Explicitly, we have 



I + V Q + iU 
Q-iU I~V 



(2.10) 



Visibilities 



Consider an interferometer with N antennas. The out- 
put signal from antenna j will be denoted £'^J^^ . As noted 
earlier, we will treat this as a two-component vector with 

components e^'^i^^, with m = X, F for a Hnear polariza- 
tion experiment or m — R, L for a circular polarization 
experiment. 

The basic datum for an interferometer is a "visibility" 
obtained by correlating a component of eout from one 
antenna with a component from another antenna: 



\ out m out n/ ' 



(2.11) 



where the angle brackets denote a time average. Both 
real and imaginary parts of this complex quantity can 
be obtained by measuring the in-phase and quadrature- 
phase correlations. 

For a fixed pair of antennas (jk), the visibilities form 
a 2 X 2 matrix V^-''^). Using equations H2.4() and H2.6() . we 
can write this matrix as 



J d^f A(^)(f)-S(r)-AW^(f) e-2"".'«-?. (2.12) 

Here the A's are the antenna patterns for the two anten- 
nas and Ujk = d'*^' — ^^^^)/A is the separation between 
the two antennas in units of wavelength. The matrix A^ 
is the hermitian conjugate of A (that is, the complex 
conjugate of the transpose of A). 

In an ideal experiment with A proportional to the 
identity matrix (no cross-polar response and identical co- 
polar response to both polarization states), V'^-''^) cx S. 
In other words, each visibility measures a simple linear 
combination of the Stokes parameters. To be explicit, let 
us define Stokes visibilities 



V, 



d^f yl(j')(f)Z(f)AW(f)e- 



(2.13) 



where Z = I, Q,U,V is a Stokes parameter. As is well 
known, these can also be written as a convolution in 
Fourier space: 



cx 



k - 27ru) 



(2.14) 



with Ajk = ^(J)^^^). 

A polarimetric interferometer can work either by in- 
terfering linear polarization states or circular polarization 
states. Throughout this paper, we will refer to these pos- 
sibilities as linear experiments and circular experiments 
respectively. Information about both linear and circular 
polarization can be obtained from either type of experi- 
ment. 

In an ideal linear experiment, we would extract the 
Stokes parameters from the visibility matrix as follows: 

Vi = ^{Vxx + Vyy), (2.15a) 

Vq = ^{Vxx-Vyy), (2.15b) 

Vu ^^{Vxy + Vyx), (2.15c) 

Vv = jiiVxY - Vyx). (2.15d) 
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Here we are assuming that all antennas split up the 
incoming radiation into orthogonal linear polarizations 
with respect to a single fixed coordinate system {X,Y). 
The superscript (jk) is suppressed. 

For the weak polarization found in CMB data, equa- 
tion (|2.15b|l is not a practical way to measure Stokes 
Q because it requires perfect cancellation of the much 
larger / contributions; in practice, such an experiment 
measures linear polarization only via U, not Q. Since 
Q U under a 45° rotation, we measure Q in practice 
by using antennas that measure linear polarization states 
in a basis (X', Y') that is rotated with respect to {X, Y). 
This can be done either by rotating the instrument or 
by having the polarizers on different antennas oriented 
in different ways. In either case, note that in general the 
Stokes parameters Q, U are not generally measured with 
the same baseline at the same time. 

The corresponding relations in a circular experiment 
are 



Vi - ^{Vrr + Vll), 
Vu = U^RL - Vlr), 



(2.16a) 
(2.16b) 
(2.16c) 
(2.16d) 



In a circular experiment, both Q and U visibilities can 
be measured simultaneously on a single baseline. 

We do not expect any cosmological source of circular 
polarization: Stokes V is expected to be zero. Nonethe- 
less, it may be useful to measure the Stokes visibility Vy 
as a monitor of systematic errors. Conversely, if a non- 
cosmological source of circular polarization is present, 
systematic errors may cause it to contribute to measure- 
ments of linear polarization. 



be represented in a circular polarization basis. The fac- 
tors R and are inserted to account for this change 
of basis. 

An ideal instrument would have J; = 1, the identity 
matrix, and Ag = A{9, 0)1. 

There is some redundancy in equation H2.17|l . Math- 
ematically, we could correctly describe any instrument 
without including 3i by simply absorbing its effects into 
Aj. However, it is convenient to maintain the distinc- 
tion between effects that happen to the signal before the 
antenna averages over the beam (effects "on the sky") 
and afterwards (effects "in the instrument" ) . Instrument 
errors are easier to model, because by definition they do 
not depend on position on the sky. 



III. EFFECT OF ERRORS ON VISIBILITIES 

In this section we compute the effects of various sorts of 
instrument and beam errors on the measured visibilities 
^Qj^u- Each of the errors considered can be modeled 
with a set of small parameters. We assume that the ex- 
perimenter has no knowledge of these errors (or else she 
would have removed them) and hence analyzes the data 
under the assumption that the experiment is error-free. 



A. Instrument Errors 

Consider first the effect of errors within the instru- 
ment, assuming for the moment that is of the ideal 
form As{r)l. We can completely characterize the instru- 
mental Jones matrix for the jth antenna with gain errors 

(i) (i) (1) (i) 

fli i52 couplings ■ 



C. Modeling systematic errors 

A wide variety of systematic errors can be modeled as 
imperfections in the matrix-valued antenna patterns of 
the antennas. We will model these errors in the following 
way: 



A(f) = J,•R•A,(f■)•R- 



(2.17) 



Here R is Rcirc for a circular experiment and is the iden- 
tity matrix for a linear experiment. 

The "instrument Jones matrix" represents errors in- 
troduced within the instrument, such as gain errors and 
cross-talk between the two outputs of a given antenna. 
The matrix As is the antenna pattern on the sky, before 
such instrumental errors are taken into account. We use 
As to model cross-polarization, beam errors, pointing er- 
rors, etc. We will always use A (with no subscript) to 
denote the overall antenna pattern including both "in- 
strument" and "sky" effects. 

We will always represent As in a Cartesian basis; that 
is, it acts on components (einx j einy)- When we are per- 
forming a circular experiment, however, ein and Cout will 



1+52 



(i) 



(3.1) 



This is similar to the characterization in ref. |17|. al- 
though our notation is not identical to theirs. In particu- 
lar, we treat the g and e parameters as complex numbers 
(to account for arbitrary phases in the errors) rather than 
introducing explicit phase angles. The parameters g and 
e will be assumed to be small (i.e., products of them will 
be neglected). 

A number of different physical effects can be encoded 
in a matrix of this form. The gain parameters g\^^ in- 
corporate both errors in the magnitude of the gain and 
unaccounted-for phase delays. The couplings ep^ can ac- 
count for mixing of the two polarization states within the 
optical and electronic systems and also, in the case of a 
linear experiment, for an error in alignment of the polar- 
izers: if the polarizers in antenna j are misaligned by an 
angle 8, then e^ '' = 5, e^^ = —5. 

By using this model for the antenna patterns in equa- 
tion (|2.12|) . we can determine the effect of all of these 
errors on the recovered Stokes visibilities. The results 
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depend on whether we are considering a hnear or circu- 
lar experiment. 



Linear experiment: Information about hnear polariza- 
tion in such an experiment comes from the visibility for 
Stokes U. Using equations H2.3|l . H2.12|l . and H2.15c|l . we 
find 



V, 



V, 



1 r 

2 



U) , A]) 



if) , Jj) 



9i 

,(fc)* 



(fc)* 
^1 

fe)* , 



(k)* 
-92 

Jk)*. 
^2 ) 



(k)* 

91 



+ 5^*) , (3.2) 



working to linear order in the small quantities. Here the 

o 

symbol V indicates the visibility that would be measured 
in the absence of systematic errors. Note that the cou- 
pling parameters e mix temperature anisotropy (/) into 
polarization; this is in general the most serious sort of er- 

o o 

ror. The Vq and Vu terms are less worrisome, since they 
only involve polarization. However, as we will see they 
do couple E to B and so can be serious for a S-mode 
experiment. 

Although for cosmological purposes we are primarily 
interested in measurements of linear polarization, with 
this experimental setup we get the visibility for Stokes 
V "for free" by subtracting rather than adding Vxy and 
Vyx [see equations (|2.15ll ]. Assuming there is no intrinsic 

o 

circular polarization {Vy ~ 0), the leading contribution 
is 



T/(ifc) _ * f/Uk) ,(3) _ 



_(]) (fe)* 



+ 4'^*), (3.3) 



neglecting terms proportional to Vq,Vu- Even very low 
levels of coupling may therefore provide a measurable sig- 
nal, which would be correlated in a known way with the 
temperature map. This may provide a useful diagnostic. 

Conversely, if there is intrinsic circular polarization 
(e.g., due to foregrounds, equation H3.2|) shows that gain 
errors can couple that signal into Vq^Vjj- 



Circular experiment: Now suppose we perform an ex- 
periment in which interference between right and left cir- 
cular polarization states is measured. Using equations 



(I?1U|) . 11321, (|2.16b|) and fTTBc|l . we find 



Uk) 



Uk) 



(fe). 



Uk)t(j) 



^Vr'{9'^ 



Vr 



+ V, 



Uk) _-^Uk) 



U 
Uk) 



(fe), 

9i 



(fe), 

92 



9i 



(fe). 



92 



u 



1 



U) 
1 



.U) 



(k)* 



-fef*) , (3.4a) 



u 



Uk), U).U).(k)* 



^2 J 



+ Vri9[^'^ 


-9f- 


(fe), 
^ .9i - 




-^V^'\9'^' 




(fe), 
-91 


+ .9f*) 






(fe)* 
— ^1 


- ^2 I 



(3.4b) 



As in the linear case, coupling errors (e) are the main 
danger, causing leakage from / into Q,U. Gain errors 
mix Vq and Vjj with each other. As we will see, this 
means that their effect on the B mode power spectrum 
is somewhat more severe than in the case of a linear ex- 
periment. 

As for linear experiments, we might consider monitor- 
ing the circular polarization visibility Vv as a diagnostic, 
even though we do not expect any cosmological signal. In 
this case, the leading term in Vy is proportional to the 
gain fiuctuations g\^^ and to V/. Since gain errors are less 
worrisome than couplings, this may not be as valuable as 
in the linear case. Furthermore, unlike a linear experi- 
ment, in a circular experiment one does not necessarily 
get Vv for free when measuring Vq, Vjj'- the linear polar- 
ization information is obtained by interfering right with 
left polarization states, while Vy comes from interfering 
identical ones [see equation H2.16|l ] . 



B. Beam Errors 

We next consider errors that can be modeled via the 
"sky" matrix A^. In this section we ignore instrument 
errors, taking to be the identity matrix. 

An ideal experiment, with equal to a scalar function 
times the identity matrix, would have identical response 
to both polarization states and no mixing between them. 
Furthermore, ideally Ag would be the same for all anten- 
nas. There are of course a large number of ways these 
idealizations can fail. Unlike instrument errors, which 
are characterized by a finite list of parameters, beam er- 
rors are characterized by arbitrary functions on the sky. 
Rather than providing a complete catalogue of this infi- 
nite space of possibilities, we focus on a few physically 
motivated possibilities. 

Beam mismatch: We first consider the case where each 
antenna pattern is proportional to the identity matrix, 
but the various antenna patterns differ from each other 
and from the form assumed by the experimenter: 



A(j)(f) = A(j)(f)l. 



(3.5) 
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In this case, we are assuming no cross-polar response and 
identical beam patterns for both polarizations in each an- 
tenna. This formulation can account for pointing errors 
as well as errors in beam shape (e.g., beam width and 
ellipticity errors). 

Assuming this form for the antenna pattern, we can use 
equation (|2.12|) to extract the Stokes visibilities, yielding 



Vr 



d?v 6-2'^^"^'= ■^Z{v)A^^^ {v)A^^> (f ). (3.6) 



For Z = {/, Q, [/, V}. These results apply to both linear 
and circular experiments, although in practice only Vu 
would be used for linear polarization information in a 
linear experiment. 

This category of error causes no leakage from / into 
polarization or even between Q and U . Nonetheless, as 
we will see it can cause E/B mixing when the power 
spectra are estimated. 



Cross-polarization: We now consider the possibility of 
cross-polar antenna response (i.e., off-diagonal entries in 
As). For simplicity, we consider only the case of an az- 
imuthally symmetric antenna. Wc therefore begin by 
determining what form of A^ are compatible with az- 
imuthal symmetry. 

First, note that A^ must be diagonal when measured 
along the x axis of our coordinate system. One way to 
see this is to invoke reflection symmetry: if ein^, cou- 
pled to eoutj, along the x axis, the coupling would have 
to change sign upon reflection, but a symmetric antenna 
should be invariant under reflections. Using polar coor- 
dinates (r, 4>) , we can therefore write 



(3.7) 

for arbitrary radial functions ^q'*, A^f' . Suppose we now 
measure the antenna pattern for some nonzero 4). We 
must get the same answer as aX (j) ~ Q, as long as we 
transform both ejn and Cout by a rotating through an an- 
gle —0 to move them to the x axis. We conclude that 
As(r, 0) = R_0 • As(r, 0) • R^, where R±0 are the ap- 
propriate rotation matrices. Performing the matrix mul- 
tiplication, we find that 



A(;) 



'4'^ + i4'^cos2</) \Afs\n2(^ ^ 
iAfsin20 -14') cos 20^ 



(3.8) 

Assuming the form H3.8|l for Ag, we obtain the follow- 



ing expressions for the visibilities: 

.i/(f )[4') {r)A^^^* (r) + A^i>^ (r)4'=)* (r)] cos 20+ 

ll/(f)[4")(r)4'=)*(r) - A['\r)A^i>{r)] sin20} , 

(3.9a) 



T/(jfe) _ 



(fv -r |{/(f )4' (r)4''^* {r)+ 



i/(f ) [4'^ {r)Af'* {r) + A'i' (r) A^"^* (r)] sin 20- 
^V{v)[A'i\r)Af>{r) - 4^")(r)4'=)*(r)] cos2</>} , 

(3.9b) 

neglecting terms that are quadratic in the small quan- 
tities Ai. Again, these results apply to both linear and 
circular experiments, although only Vjj is used for po- 
larization measurements in a linear experiment. We im- 
mediately see that cross-polarization has the danger of 
coupling / to polarization. 

If we write down a similar expression for the circu- 
lar polarization visibility Vy , we find terms proportional 
to Stokes Q and U but not /; thus no great insight into 
cross-polarization is likely to be found by monitoring Vy ■ 
On the other hand, if there is a strong intrinsic circular 
polarization signal from foregrounds, cross-polarization 
may cause it to contaminate the linear polarization ob- 
servables. 



IV. EFFECTS ON POWER SPECTRA 
A. Introduction 

The primary goal of almost any CMB experiment is 
to measure some or all of the temperature and polariza- 
tion power spectra. We must therefore consider how to 
propagate the errors described above to obtain forecasts 
of the effects of various systematic errors on power spec- 
trum estimates. The completely correct approach to this 
question is to define a precise experimental setup and 
simulate it in detail. We will instead adopt an approach 
that is simpler and more general. 

The power spectra are well-localized in Fourier or 
spherical harmonic space and therefore not at all local- 
ized in image space. Interferometric visibilities are fairly 
well-localized in Fourier space, so power spectrum esti- 
mates can be obtained from individual visibilities (or in 
the case of polarization, from pairs of visibilities Vq, Vu). 
This is in contrast to an imaging experiment, in which 
the entire map goes into each power spectrum estimate. 
We can therefore assess the effects of various errors on 
power spectrum estimates by working with one visibility 
pair at a time. Naturally, analysis of real data would 
be considerably more sophisticated than the process we 
consider here, but this approach is adequate for assess- 
ing the rough scale at which each systematic error affects 
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power spectrum estimates. The rest of this section will 
present this approach in detail. 

A data set from an interferometric experiment consists 
of measurements of the visibilities Vq, Vjj for many dif- 
ferent baseline separations and instrument orientations 
- that is, for many different values of u. Let us as- 
sume that our data set consists of a set of visibilities 
Vqi , Vc/i , Vq2 , Vu2 Vqn , VuN ■ We will assume that 
Vqi, are both measured with the same baseline vector 
Ui, and furthermore that distinct baselines u^, are far 
enough apart in the Fourier plane that we can treat the 
corresponding visibilities as uncorrelated: ju— Uj| > Au, 
where the Fourier-space resolution Am scales inversely 
with the beam width. In a real experiment, these as- 
sumptions would presumably not be true: at least some 
regions of the visibility plane would be oversampled, and 
in the case of a linear experiment the two Stokes pa- 
rameters would not always be measured with identical 
baselines. In such a case, we can imagine binning the 
visibilities so as to consider a smaller number of effec- 
tively independent samples in the visibility plane. 

In this approximation, each set of visibilities (Vgi, Vjji) 
gives us an independent estimate of the power spectra 
Cf",Cf at I « 2TTUi. (Here "independent" means "inde- 
pendent of the estimates obtained from the other visibil- 
ities" ; unfortunately, it is not in general the case that the 
estimates of and are independent of each other.) 
We therefore begin by assessing the effect of each system- 
atic error on the power spectrum estimates derived from 
a single pair of visibilities {Vq, Vjj) for some fixed base- 
line vector u. We assume that systematic errors must be 
controlled at least well enough that the fractional errors 
introduced in the power spectrum estimates from each 
individual visibility are small. 

In this section, we begin with the optimal estimators 
of E and B band powers from a single visibility pair 
{Vq,Vu) under the idealized assumption that there are 
no errors in the experiment. We then imagine "turning 
on" sources of systematic error, one at a time, and show 
how to calculate the resulting root-mean-square (r.m.s.) 
band power errors. We will find that the r.m.s. error 
induced by any particular systematic error can be written 
in the form 

{SC^^^f = pI,, kI.^^C'c'. (4.1) 
I J 

Here p is a parameter characterizing the strength of the 
source of error, (e.g., one of the gain fluctuation param- 
eters gf^)- The superscript K — {E,B} indicates the 
type of power spectrum being measured, and ^C^^g is 
the error in an estimate of a band power. The quan- 
tities C^,C'^ are band powers with I, J ranging over 
{T, X, E, B} (temperature, TE cross-correlation, and E 
and B polarization spectra) . We will see how to calculate 
the coefficients k below. 

As Figure ^ shows, there is a clear hierarchy in the 
input power spectra: > C'^ > > . The above 
sum is therefore usually dominated by one or at most two 



terms that couple large spectra to smaller ones. For each 
of the sources of error described in the previous sections, 
we can find the one or two most serious such couplings 
and thus estimate the level to which the overall error 
must be controlled or at least understood. 

As noted above, our approach is to demand that errors 
be understood to a level low enough that they do not 
contaminate band powers measured from each individ- 
ual baseline. This may be regarded as a conservative ap- 
proach: even if an error has quite a large effect on power 
spectrum estimates from a single baseline, it may still 
be possible to extract good power spectrum estimates by 
combining results from a set of many baselines. This is 
particularly true if the visibility plane is heavily oversam- 
pled, with many different visibilities within a resolution 
element Am, if the errors associated with the various vis- 
ibilities are independent and average to zero. However, 
since systematic errors are in general non- Gaussian and 
time- varying, one would not want to rely on such an av- 
eraging procedure unless the statistical properties of the 
errors are understood quite well. It therefore seems rea- 
sonable to adopt such a conservative approach. 



B. Ideal estimators 

We begin with the estimators for E and B band pow- 
ers for a single visibility pair {Vq,Vu), for an idealized 
experiment with all systematic errors "turned off." In 
the following subsection, we will examine the effect on 
these estimators when the errors are introduced. 

Let V = {Vi, Vq, Vu) be a set of visibilities correspond- 
ing to a single baseline u. Our hypothetical experiment 
measures Vq and Vjj, but not necessarily V/. As we will 
see, it is convenient to include Vi in the formalism any- 
way, because it may be coupled to the others by various 
systematic errors. 

In the absence of systematic errors, the / visibility 
measures only the temperature fiuctuations: 

{\Vi\')=Cl^. (4.2) 

Here is a band power giving the average tempera- 

ture power spectrum in the vicinity of I = 2ttu. We omit 
a constant of proportionality in this expression and here- 
inafter by assuming that all visibilities have been scaled 
by an appropriate factor. 

Suppose for the moment that the baseline vector u 
points in the x direction of our chosen coordinate system. 
The contributions of E and B polarization to Stokes Q 
and U are expressed in Fourier space as follows: 

g = i?cos(2(/)) - Ssin(2(/)), (4.3a) 
[7 = i?sin(20) -|-i?cos(20), (4.3b) 

where is the angle made by the wave vector k with 
respect to the a; axis. Each visibility is an average of Q 
or U over a range of wave vectors centered on k = 27ru 
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FIG. 2: The quantity as a function of antenna separation. 
The quantity on the horizontal axis is the product of baseline 
and beam width and is proportional to the separation be- 
tween the antennas. A Gaussian beam shape is assumed. For 
antennas of diameter D with FWHM= 1.22A/D, the vertical 
dotted line at ua = 0.52 corresponds to a pair of antennas 
that are just touching. The dashed curve is the approxima- 
tion il^Tli . 



[see eq. (|2.14ll ]. Assuming that the power spectra can be 
pulled out of the integral to give band powers, we have 



{\vuf) = a 



(4.4a) 
(4.4b) 



where , are averages of cos^ (2</)) , sin^ (20) over the 
antenna patterns: 



^_ /|yl2(k-27ru)|2cos2(20)d2k 
/|22(k-27ru)|2d2k 

— _ /|l2(k-2™)|2sin2(20)d2k 



/|A2(k)|2d2k 



= 1 



(4.5a) 
(4.5b) 



Here is the Fourier transform of the squared antenna 
pattern. There can also be T-E correlations, given by 



The quantity characterizes mixing of E and B 
modes within each visibility pair. As a result, the de- 
gree to which systematic errors couple different power 
spectra is strongly dependent on . Fig. |21 shows 
as a function of ucr, the product of baseline length and 
beam width. This product is also proportional to the 
antenna separation in units of antenna diameter. In par- 
ticular, for antennas with a Gaussian beam pattern with 
FWHM = 1.22A/-D, the separation is equal to the diam- 
eter (the antennas are just touching) when ua — 0.52. 
The quantity is well approximated by 

(4.7) 



27r2(ua)2' 

Given the visibilities Vq, Vjj , the optimal estimators of 
the polarization band powers are 



(7^=7(c2|l^Q|2-s2|yc/n, 

C^=7(?I^4/|'-?|Fq|2), 
suppressing the subscripts 27rM, with 



(4.8a) 
(4.8b) 

(4.9) 



the cross-correlation . These relate Vq to V/: 



One can check using the covariances above that these 
expressions give unbiased estimates. Furthermore, it is 
straightforward but tedious to check that C^, are the 
maximum-likelihood estimators in the case of Gaussian 
fluctuations. The Cramer- Rao inequality (e.g., jl^l) then 
implies that they are the optimal estimators. 

It will be convenient to write these estimators in fancier 
linear-algebra language. The covariances of the Stokes 
visibilities v = (V/ , Vq, Vj/) can be organized into a 3 x 3 
covariance matrix: 

(C^ _C^c _ 
M, = (vvt) = C^c C^c2 + C^s2 _0 _ 
V C^s'^ + C^c'^j 

(4.10) 

The estimators of and are quadratic functions of 
the vector v: 

C'^ = vt.NE-v, • Ns • V, (4.11) 



r 



(4.6) 



where c is an average of cos(2(/)) analogous to equations 
(|4.5|) . These expressions are valid only for u parallel 
to the X axis. Without this assumption, there would 
be a Vj-Vu covariance proportional to s, the average of 
sin(2(/)), and a Vq-Vu covariance proportional to sc, but 
with this geometry both of these terms vanish. 

In the limit where the visibility is measured with two 
very widely separated antennas {2ttu(t ^ 1 for beam 
width a), each visibility samples a very narrow region in 
the Fourier plane. In this case, c2 « 1, c ~ 1, and s"^ ~ 0. 
The squares of Vi, Vq, Vu then provide pure estimates of 
C'^,C^,C^ respectively. 



where 

/O 0_ \ /O 0_ 0^ 

Ne=7 0c2 0_ , Ns=7 0-^0_ 
\0 -s^J \0 c2^ 

(4.12) 

In this notation the proof that is unbiased looks like 
this: 

(C^) = (vt . • v) (4.13a) 
= Tr {Ne ■ (vvt)) (4.13b) 
= Tr(NB-M„) (4.13c) 
= C^, (4.13d) 
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using the explicit forms of M^, and N^; in the last step. 

Throughout this section, we have made the simplifying 
assumption that the baseline u was parallel to the x axis. 
We now generalize the results to the case where u points 
in an arbitrary direction. Let a be the angle between 
u and the x axis. We can calculate the optimal power 
spectrum estimators by rotating to a coordinate system 
in which u is on the axis before applying the above pre- 
scription. The rotated Stokes vector is 

/I \ 

Vrot=R-v= cos2q! sin2a -v. (4.14) 
\ — sin 2a cos 2a J 

The covariance matrix (|4.1UI) applies to the rotated vec- 
tor Vrot ; the original unrotated vector thus has covariance 
matrix (v • vt) = R-^ (v^t • vj^^) • R = R^^ • M„ • R. Sim- 
ilarly, when applied to the unrotated data the matrices 
'Ne,b are simply replaced by R^^ • Ne,b ■ R- 

In the next section we will continue to examine the spe- 
cial case of u pointing in the x direction for simplicity, 
but these transformations can always be made to gener- 
alize the results. 



C. Introduction of errors 

Suppose that there is an error (e.g., a gain error) in 
the data, but we don't know it. Since we think we are 
dealing with an error-free experiment, we use the opti- 
mal prescription (|4.8I) to estimate the power spectra from 
a single set of visibilities. The presence of the system- 
atic error will alter the covariances (|4.10() and hence the 
statistical properties of the estimators. Let 6C^ with 
K = {E, B} be the difference between the estimate we 
actually get and what we would have gotten in the ab- 
sence of systematic errors: 

SC — C'actual ~ Aio errors- (4-15) 

It is natural to use the r.m.s. values of these differences, 

<5A^„3 = {{SC^ry/^ (4.16) 

to quantify the effect of each systematic error on the 
power spectrum estimates. We focus for the moment on 
the case of "instrument errors" as described in Section 
nil Al deferring the generalization to beam errors until 
section IVll 

The effect of each instrument error is to mix together 
the Stokes visibilities Vj,Vq, Vjj in a linear fashion (ne- 
glecting circular polarization for the present). If v de- 
notes the vector of Stokes visibilities that would be ob- 
tained without systematic errors, the vector actually ob- 
served is 

v'=v + (5v = v-|-E-v (4-17) 
for some 3x3 matrix E. 



The error in the power spectrum estimates is 

SC^ = (v + (5v)t • Ni^ • (v + (5v) - vt • Na' • v. (4.18) 

Assuming that the errors are small, we can neglect the 
term that is quadratic in 5v. 

SC^ = • (Et • + Na' • E) • V = vl' • • V. (4.19) 

Assuming Gaussian fluctuations, there is a relatively sim- 
ple expression for the variance of this quantity, as shown 
in the Appendix: 

((54^3)2 = Tr[(AA • M,)2] + [Tr(AA • M,)]\ (4.20) 

As noted in the previous subsection, when these for- 
mulae are applied to visibilities in a coordinate system 
that is not aligned with the baseline vector, we must cor- 
rect them by conjugating with the matrix R. All that is 
necessary is to replace the error matrix E with 

E^R-E-R-^ (4.21) 

For any particular source of error, we now have a recipe 
for calculating the effect on the power spectrum: we write 
down an explicit form for the matrix E and apply equa- 
tion (|4.20|1 . The result will be a sum of terms that are 
quadratic in the band powers C'^ . All of 

these quantities represent band powers at the same mul- 
tipole I = 27ru; we continue to omit the multipole sub- 
script 2t:u throughout this section for simplicity. In par- 
ticular, the error in our estimate of the i?-mode band 
power generically looks like 

- V^AC^f + VrxC^C" + v^eC^C + ... 

+ rj^EiC^r + .... (4.22b) 

Each of the coefficients 77 depends on the various param- 
eters that characterize the instrument errors, such as the 
gain fluctuations gp'' and couplings ep-*. Of course, a 
similar expression would apply to . 

We now imagine "turning on" one error at a time. Con- 
sider a systematic error characterized by a single param- 
eter p. The coefficients ry in the above expression will 
contain terms proportional to at leading order, be- 
cause equation (|4.2Uf) is quadratic in the error matrix E: 

Vu = <.,.,,pLs- (4.23) 

Here K = {E, B} is the power spectrum we are trying to 
estimate; and I,J = {T, X, E, B} are the power spectra 
being coupled to our estimate. 

Equations 14.22b|l and H4.23|l together yield the key 
result of this section: 

I, J 
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FIG. 3: The coefficients 



for iinear-experiment gain errors, plotted as a function of polarizer angle a and antenna 



separation ua. The left plot is for the parameter Re(72) 
and so is not plotted here (but see Fig. 0J. 



and the right is for 101(72). The coefficient for 71 has no Q-dependence 



We will drop the subscript rms on p below. 

For any given type of error, of course, some of the 
coefficients 77 will vanish. For instance, as we saw ear- 
lier, gain errors do not couple / to Q, J7, so there will be 
no contributions coupling or to the polarization 
band power estimates. Since there is a clear hierarchy 
> > > (see Fig.lH), it often makes sense 
for each error to consider only the term in 14.11) that con- 
tains the biggest power spectra. 

For any one source of error, the fractional error on the 
band power will be 



K 
rms 



VcFcJ 



(4.24) 



assuming that one term in the sum (|4.1|) dominates the 
error. If we demand that this fractional error be below 
some specified tolerance, then we can determine the re- 
quired specification for the input parameter p. The coef- 
ficients K are thus the key to assessing the severity of any 
particular source of systematic error. The next sections 
present calculations of these coefficients. 



V. RESULTS: INSTRUMENT ERRORS 

This section presents the results of applying the above 
formalism to the various instrument errors described in 
section lTlI Al Beam errors will be treated in the following 
section. There are several cases to consider. A table 
summarizing the key results, along with a discussion of 
their implications, may be found in Section [VIII 

Gain errors: Linear experiment: Consider an exper- 
iment that measures linear polarization states, and as- 
sume that there are gain errors gp"* , ignoring couplings 

(ep-*) for the present. As equation (|3.2|l shows, the re- 
sulting error in each visibility is simply proportional to 
the visibility itself: SVq = jqVq and SVu = juVu- Here 



1/ 0") , (j) , (fc)* , (fe)*N 

1Q,U = 5(51 + 92 +9l +92 ) 



(fe)* 
92 



(5.1) 



with the parameters gp-* evaluated at the time the corre- 
sponding visibility is measured. If Vq,Vu are measured 
with the same antennas (by rotating the polarizers in 
each antenna 45°), and if they are measured at nearly 
the same time so that the gains have not drifted, then 
7q = lu ■ It is far more likely, however, that the two vis- 
ibilities have independent gain fluctuations, in which case 
7Q and 7(7 should be treated as independent, unknown 
error parameters. 

In the notation of the previous subsection, wc can char- 
acterize these errors with a matrix 



^0 
E = I 71 + i72 

.0 71 - 



where 



71 = \{lQ+lu)-, 
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lu) 



(5.2) 



(5.3) 



In general, we will concern ourselves only with errors 
that couple larger power spectra to smaller ones. In an 
experiment to measure E modes, gain fluctuations do not 
lead to any such terms. We therefore focus on a i? mode 
experiment. If we use this matrix to calculate the errors 
on the B power spectrum, the leading term is the EE 
term: 



E\2 



(5.4) 



The coefficient is 



■[|72|^sin2(4a)-|- 



2(c2 - s2)2 

4i2^(4Re(7i)2 + (3-4sin2(4a))Re(72)^)]. (5.5) 
The baseline vector u makes an angle a with the x axis. 
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FIG. 4: Coefficients k averaged over a for linear-experiment gain fluctuations. The left plot shows for the three 

independent gain fluctuation parameters Re(7i), Re(72), Im(72) as defined in equation 15.311 . The right plot shows the coefficients 
Kg . In each case, the dashed curves show the leading term in a Taylor series in s^, which is generally a good approximation. 
The vertical dotted line at ua — 0.52 corresponds to antennas that are touching. 



The effect is characterized by three parameters jir = 
R-e(7i), 72r = Re(72), 72, = Im(72). The coefficients 
associated with these parameters are 



(s2 c2)2 



(c2 ~ s2)2 ■ 



(5.6a) 



kI = — ^— [sin^(4a)- 



4(3cos2(4a) - sin2(4a))s2 c^], (5.6b) 
s2 c2 sin^(4a) 



2(c2 - s2)2 



(5.6c) 



We can simphfy these expressions and those to follow 
in two ways. First, since a typical experiment will in- 
volve visibilities measured with many different baseline 
orientations, we will generally average over the angle a. 
Second, since s'^ is generally a small quantity, we can of- 
ten keep only the leading term in a Taylor expansion in 
s2. In these approximations, the coefficients simplify to 



'^B.EE.-ilr ~ 



(5.7a) 
(5.7b) 



Figure |3 shows the coefficients associated with 72 as 
functions of both a and antenna separation ua (which 
determines s2). Figure 01 shows the result of averaging 
over a, as well as the 71^ coefficient, which is indepen- 
dent of a. As the figure indicates, the leading-order ap- 
proximation in s2 ig quite good. 

We should consider whether it is always adequate to 
keep only the leading term rj^^- The full expression for 
the power spectrum error contains a term rj^g as well. 
T he fracti onal error caused by this term scales as only 
y/C^JC^ rather than C^/C^ (ea. KWi . However, the 
coefficient can be small (see Fig.^J, especially for 

widely separated antennas. As the right panel of Fig. ^ 



shows, the coefficient kb.eb can be much larger than 
hb,ee- Determination of which term dominates must 
unfortunately be made on a case- by-case basis. 

The size of the coefficient k is largely determined by 
its dependence on the small parameter . As a general 
rule, we need to weigh the importance of a subdominant 
contribution to the error (i.e., one that ranks lower in 
the T, AT, E, B hierarchy) if that error has a weaker de- 
pendence on s2. In the case of the parameter 71, for 
example, the EE term has oc but oc \/~^ . 

For 72, the EB coefficient has an s2_iudependent term. 
Thus for small (large separation), the EB terms may 
become more important than the EE terms. See Section 
IVIII and especially Table U for further discussion of this 
point. 

Gain errors: Circular experiment: In this case, we as- 
sume an experiment that measures Vq and Vjj simultane- 
ously by interfering right and left circular polarizations. 
Again, we focus on an experiment to measure B modes, 
as gain errors do not pose a serious problem in an E mode 
experiment. From equations H3.4|l . we see that gain errors 
f;P^ produce an error matrix of the form 



where 



^0 
E = I 71 «72 I , 

,0 -^72 7i 



1 1 (i) I (i) I (fe)* I 

71 = 5(51+32 +51 +52 )> 

72 = 5(51 52 -ffi +52 )■ 



(5.8) 



(5.9a) 
(5.9b) 



As in the previous case, the dominant error contribution 
to a measurement of B power is the EE term: 



B _ 2s2c2(|72|2+4s2c2(Re(7i)2-Re(72)')) 
Vee — 



(c2 - s2)2 



(5.10) 
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FIG. 5: Coefficients for gain errors in a circuiar-poiarization experiment with parameters defined in (15.911 . As in Fig. |3] tfie 
ieft piot sfiows tlie coefficients and tfie riglit shows Dashed lines indicate the leading-order approximation in s^. 

These coefficients are independent of a, so no averaging was necessary. The vertical dotted line corresponds to antennas that 
are touching. 



The coefficients are 



(s2 c2)2 



2^2 



(c2 - s2)2 



(^2) 



2s2c2(l-4s2c2) „^ 



(c2 - s2)2 
2^2^ 



(c2 - s2)2 



2s2 



(5.11a) 
(5.11b) 
(5.11c) 



where the approximate equahties are the leading terms 
in an expansion in s2. See Figure |5| 

As in the case of a linear polarization experiment, 
the EB error term can become dominant for widely- 
separated antennas. In particular, for the parameter 
72 the EB coefficient has a term independent of 52^ 
'^B EB ^ ^ leading order in s'^. 

Couplings: Next we turn to errors parameterized by 
the "coupling" terms e" in the instrument Jones ma- 
trix. These errors include electronic cross-talk as well as 
errors in the alignments of the polarizers in a linear ex- 
periment. In both linear and circular experiments, these 
errors couple / into Q, U, so the dominant terms will be 
those involving the temperature power spectrum. In this 
section, unlike the previous ones, we should consider E 
as well as B mode experiments. 

The error matrix characterizing I Q,U leakage in 
this situation is 



E 




(5.12) 



Here £1,62 are the coefficients of Vi in equations (|3.2() 
and H3.4|) : For a linear experiment. 



£1,2 



(4 



(5.13) 



ii) 

with the parameters el evaluated when the correspond- 
ing visibility is measured. For a circular experiment. 



£1 - al^i 



£2 



(k)* 



(5.14a) 
(5.14b) 



In the linear case, there are also terms that couple Q and 
U . We omit these, as the errors they produce are always 
small in comparison to the terms involving /. 

Consider first a B mode experiment. The TT and TX 
terms in (j4.22b|l vanish, so the dominant contribution is 
the TE term. After averaging over a, this term is 



B 



'C2(|ei|2 + |£2|2) 



(c2 - s2)2 

SO the coefficients for 
Re(ei), Im(ei), Re(e2), R.e(e2) are 



the 



=2^2 



(5.15) 
parameters 

(5.16) 



as shown in Figure El 
independent piece: 



The TB term has an s^- 



_ 1 + 3s2 c2 
^ (^-^)2 ~ ' 



(5.17) 



which can be important for large antenna separation. See 
Table n 

The coefficients for an E mode experiment are the 
same as for a B mode experiment with E and B switched. 
The dominant contribution is therefore « 1. 



VI. BEAM ERRORS 

In section lTV CI we derived a method of forecasting the 
effects of instrument errors on power spectrum estimates. 
We now generalize this method to the case of beam errors. 
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FIG. 6; Coefficients for coupling errors parameterized as in 
equations 15.1311 and 15.1411 . Dashed and dotted curves are as 
in the previous figures. 



In the case of instrument errors, we were able to write 
the errors in the visibihties as (5v = E • v, where the error 
matrix E depended only on the unknown error parame- 
ters. Beam errors cannot be treated in this way. Both Sv 
and V are integrals over the Stokes parameters, but with 
different weightings in Fourier space. As a result, one is 
not a simple linear transformation of the other. Fortu- 
nately, for a number of important sources of error, the 
differences in Fourier-space weighting are modest: the 
errors sample roughly if not exactly the same regions of 
the Fourier plane. We can therefore still express our fi- 
nal results in the form H4.1|l . after we have made some 
adaptation to the formalism of section HV CI 

Combine v and (5v together into a 6-dimensional vector 
w = (v,(5v) = {Vi,Vq,Vu,SVi,SVq,SVu). To leading 
order in (5v, the error in the power spectrum estimate is 



6C 



K 



Mi 



K ■ W, 



(6.1) 



where K = {i?, B} and the matrix M can be written in 
block form as 



Nk 



N 



K 







(6.2) 



with N/f the same as for instrument errors. It is straight- 
forward to check that this reduces to H4.18|l . 

Using the identity proved in the Appendix again, we 
can write 



K \2 
rms/ 



Tr [{Nk ■ M^)2] + [Yr{NK ■ M„,)]' , (6.3) 



where M,„ = 



r^) is the covariance matrix of the 



vector w. We now need a recipe for calculating the ele- 
ments of this covariance matrix, which will contain terms 
proportional to the various power spectra. 



In sections III Bl and IIII Bl we expressed each compo- 
nent of V and 5-v as an integral over the Stokes param- 
eters. To be explicit, let s(k) = (/(k), (3(k), t/(k)) be a 
vector giving the Fourier-space Stokes parameters. Each 



component of w can be expressed in the form 

j d2fcW,(k) •s(k) (6.4) 

for some vector- valued window function W^. A covari- 
ance matrix (wiW*) then becomes an integral over k of 
the two window functions times the covariances of the 
Stokes parameters (s • s'''). The latter are proportional to 
the input power spectra. So once we have written down 
the window functions for the visibilities and their asso- 
ciated systematic errors, we can calculate an expression 
giving contributions to the error SC^^ in terms of the 
input power spectra, just as in the case of instrument 
errors. For some parameterizations of beam errors, the 
resulting integrals can be performed analytically to yield 
closed-form expressions like those in the previous section 
for the coefficients k. However, the resulting expressions 
are complicated and unenlightening, so we present the 
results of numerical integration instead. 

As in the previous section, we now examine detailed 
case-by-case results. Section IVIII provides a summary of 
the implications. 

As noted earlier, the set of possible forms for beam er- 
rors is dauntingly large. Our treatment will necessarily 
be restricted to a small set of physically motivated possi- 
bilities rather than exploring the entire space. As in the 
previous section, we will imagine turning on one error at 
a time. 

Differential pointing errors ("squint"): Suppose that 
some antennas have slight pointing offsets relative to oth- 
ers. This situation can be treated as a beam mismatch 
error as in equation H3.5II . 

Let Ao{r) be the antenna pattern in the absence of 
the pointing errors, which we will take to be a Gaus- 
sian. Then in the notation of equation H3.5|l . the antenna 
pattern for the jth antenna is 



A'^J^r) = A(f + <5fj). 



(6.5) 



where drj is the pointing error of the jth antenna. Using 
equations (|3.6|) . we find that each visibility looks like 



Vz = 



re 



Z(r)A(^'(f)A(^-)*(r), (6.6) 



where Z = {Q,U}. The product of two Gaussians 
is a Gaussian centered at the midpoint of the two: 
^b)(f.)yl(fc)(f.) oc exp[-(f- i((5fj- +(5ffc)2)/(2CT2)]. That 
is, each visibility is calculated using an effective beam 
pattern that is shifted by the average of the shifts of the 
two antennas. For any given antenna pair (jfc), we define 
an error parameter 



2a '■ 



(6.7) 



the average of the two antennas' pointing errors in units 
of the beam width. 
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FIG. 7: Error coefficients for beam pointing errors. The quan- 
tity 6+ is a common pointing error (the same for both Vq,Vu), 
while S- is a relative pointing error. Both S± are measured in 
units of the beam width. Resuhs are averaged over directions 
of S±. Furthermore, in the case of S-, an average has been 
taken over the angle a between u and the x axis. Dashed and 
dotted curves are as in the previous figures. 

Shifting a function by an amount A is equivalent to 
multiplying its Fourier transform by e*'' '^, so by (|2.14|) 
the visibility is 

Vz^ J dPk ^(k)l2(k - 27ru)*e"'(''-2^") ''^'''". (6.8) 

To leading order in 8jk^ the error is 

5Vz = ~i j d^A:Z(k)^(k-27ru)*(k-27ru)-5jfecr- (6.9) 

As before, we imagine a single measurement pair 
(Vq,Vc/) corresponding to the same baseline u. Let 5q 
be the value of 5jk corresponding to the visibility Vq and 
5u be the value corresponding to Vjj. In a circular ex- 
periment, the two visibilities are measured with the same 
antenna pair, so 5q = du, while in a linear experiment 
they should be regarded as independent error parame- 
ters. It is convenient to express our final results in terms 
of the sum and difference, 

d± = ^iSQ±Su). (6.10) 

In a circular experiment, ^_ = 0. 

We can now calculate the various correlations 
(SVqVq), etc., by integrating over k. The result will con- 
tain terms proportional to the band powers C^„, C^„, 
and quadratic in the parameters (5+ , <5_ . We can there- 
fore define parameters k characterizing these errors ex- 
actly as in the case of instrument errors [equation (|4.1(l ]. 

As in the previous cases, the severity of the errors de- 
pends on , which characterizes the degree of EB mixing 
within each visibility pair. Furthermore, the components 
of S± parallel and perpendicular to the baseline, Sj.^^ and 
6±±, contribute differently. Finally, in the case of 
the results depend on the angle a between the baseline 



and the coordinate axis. For simplicity, we have averaged 
over both components and d± for each of d±, and we 
have also averaged over a. Figure illustrates the re- 
sulting coefficients. As before, the coefficients are well 
approximated by the leading-order terms in an expan- 
sion in s^, which are given in Tabled Not surprisingly, 
J_ is a greater source of error than d+. As comparison 
with Figures 21 and O indicates, the effects of differential 
pointing errors (d-) are generally similar to those of gain 
errors. 

Beam shape errors: Equation (|3.5|l can also be used to 
model errors in the beam shape. To illustrate this, we 
consider Gaussian beams with errors in the beam width. 

Assume that in an ideal, error-free experiment all an- 
tennas have azimuthally symmetric Gaussian beam pat- 
terns with beam width a. Suppose that in actuality 
each antenna has an elliptical beam pattern with differ- 
ent beam widths ctP\(T2"''' along its two principal axes. 
As equation H3.6I) indicates, the effective beam pattern 
for each visibility F^-''^) is just the product of the two an- 
tenna patterns. The product of Gaussians is a Gaussian, 
so the effective visibility beam pattern will be of the form 

^(j)(f.)^(fc) (f) ^ exp[-f • (1 + AjTc) • f/(2cT2)]. (6.11) 

Here the symmetric 2x2 matrix A characterizes the de- 
viation of the beam from the ideal symmetric Gaussian of 
width a. To be specific, the eigenvectors of Ajk give the 
two principal axes of the elliptical beam, and the beam 
widths are o j \J\ + \i where Ai,A2 are the eigenvalues. 
We will assume that the errors are small and work to lead- 
ing order in A^. The fractional errors in the beam width 
in the two principal directions are then buija — —Xijl, 
where ? = 1, 2 label the two principal axes of the beam. 

As usual we consider a visibility pair (Vq,V;7) mea- 
sured with a common baseline u. In the case of a lin- 
ear experiment, the two visibilities may be measured 
with different antenna pairs, so we should consider two 
sets of beam shape parameters characterized by matrices 
Aq, A^/. In a circular experiment where both visibilities 
are measured simultaneously with a single antenna pair, 
Aq — Aj/. As we have seen before, we can treat both 
cases simultaneously by defining 

A± = i(AQ± Ac/). (6.12) 

The matrix A_|_ characterizes the average beam shape 
when the two visibilities are measured, and A_ charac- 
terizes errors in beam shape that differ between Vg , . 
In both cases, the two eigenvalues of the matrices give 
fractional errors in beam width in the two principal di- 
rections: 

C±,* = = ^- (6.13) 

Here i = 1,2 labels the two principal axes for each of 
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FIG. 8: Error coefficients for beam shape errors, parameter- 
ized by the coefficients ("common") and ("difF.") in 
equation (Iti.l3ll . Dashed and dotted curves are as in previous 
figures. 



We will refer to errors parameterized by C+^i as com- 
mon beam shape errors and to those parameterized by 
as differential errors. For a circular experiment we 
expect differential errors to vanish, while for a linear ex- 
periment both should be of comparable magnitude. We 
will consider separately the effects of common and differ- 
ential errors. 

In each of the two cases, there are three parameters: 
Ci,C2j and an angle /3 giving the orientation of the prin- 
cipal axes relative to the coordinate axes used to define 
Q, U. For common errors, the results are independent 
of /3, but for differential errors they depend on /3 (unless 
= C-,2, in which case there is rotational symmetry). 
We assume that the principal axes are randomly oriented, 
so we average over f3 in the results below. 

For both common and differential errors, there are two 
qualitatively different possibilities one might wish to con- 
sider. If Ci = ^2, then the beam is circular, and we have 
made an error only in its width. On the other hand, 
the case Ci = — C2 corresponds to a pure beam shape er- 
ror, with the beam stretched along one axis and squeezed 
equally along the other. Of course, the most general case 
would be a combination of the two. The final results (af- 
ter averaging over (3 where appropriate) turn out to be 
the same in both cases: the error depends only on the 
combination C,f + Cl regardless of the relative signs. 

Figure |S1 shows the coefficients k associated with beam 
shape errors. The results are quite similar to those for 
pointing errors. In particular, the differential errors that 
arise in a linear experiment are more severe than the 
common errors, which arise in both linear and circular 
experiments. 

Cross-polarization: The final case we consider is 
azimuthally-symmetric cross-polarization, with antenna 
patterns of the form (|3.8|) . We consider an ideal exper- 
iment to be one with cross-polar terms A^^^'' = for all 
antennas. The error term can in principle be an arbitrary 



FIG. 9: Coefficients for cross-polar beam response, parame- 
terized by jj,Q , [eq. 16.15ll l . Dashed and dotted curves are 
as in previous figures. 

function of r. We generally expect cross-polar response 
to be small near the beam center, so we adopt the fol- 
lowing simple form for the cross-polar response: 

4\r) = ^^.,^Ao{r), (6.14) 

where Aq is assumed to have the usual Gaussian form and 
fj,i is the parameter characterizing the size of the error. 

As an aside, note that this particular form arises in one 
simple model of an antenna. Suppose the antenna lies in 
the xy plane and responds equally to both x and y com- 
ponents of the incoming electric field, with no sensitivity 
to the z component. In the fiat-sky limit such an antenna 
has no cross-polar response, but sky curvature introduces 
cross-polarization of this form (because Eg is reduced a 
factor of cos9 w 1 — ^9^ upon projection onto the xy 
plane, while is unchanged). This cross-polarization 
is characterized by /i = icr^ with a in radians. (Inci- 
dentally, when sky curvature is taken into account one 
must be careful to distinguish among inequivalent defi- 
nitions of "cross-polarization." The most natural one in 
this context, because it respects azimuthal symmetry, is 
"definition 3" in ref. 

The relevant quantity for characterizing the error in 
each of the visibilities Vq , Vu is 

MQ,(7 = (6-15) 

the average of the two fi parameters when each of Vq , Vu 
is measured. As usual, for a circular experiment /ig = fiu 
while in a linear experiment the two are independent. In 
this case, however, it makes no difference which case we 
consider, as the error contributions due to /iq, simply 
add independently (in quadrature). 

Figure |51 shows the leading error coefficients for this 
case. Since these errors couple / into Q, U, the dominant 
terms are those involving the temperature power spec- 
trum, and errors can be quite significant for both E and 
B measurements. 
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Experiment Measurement 
type 


Error 
source 


Dominant 
contribution 


Fiducial 

sc/c 


Secondary 
contribution 


Fiducial 


Linear 


B 


Gain error*^ 




2 


2172 




_ 1 

^ 2 


8.772 


Circular 


B 


Gain error'' 




\/2p 


6072 




= V2 


2472 


Linear/Circular 


B 


Goupling^ 






730£ 


^B,TB,e 


= 1 


300e 


Linear/Circular 


E 


Coupling^ 




= 1 


17e 








Linear 


B 


Pointing'^ 


'^B.EE,5_ ^ 




30(5_ 


'^B,EB.S_ 


= 1 


17(5_ 


Circular 


B 




'^B,EE.Sj_ 






'^B,EB.Sj_ ^ 






Linear 


B 


Beam shape 


'^B.EEX- ^ 


\/p 


42C- 


^B.EBX- 


= 1 


17C_ 


Circular 


B 


Beam shape 


^B,EEX + 


3.5s2 


21C+ 


'^B,EBX+ ^ 


Vss^ 


7C+ 


Linear/Circular 


B 


Cross-polarization^ 


^B.TE.^L 


\/P 


730/i 


^B,TB,fi 


= 1 


300^ 


Linear/Circular 


E 


Cross-polarization^ 


'^E.TE.fi 


= 1 


17Ai 









"Eq. EHJ. 

'^Eq. I5 J3I (linear); Eq. 15.141 (circular). 

TABLE I: Effects of instrument errors (above line) and beam errors (below line). See Section fVllI for details. 



VII. DISCUSSION 

This paper has presented a method of quantifying the 
effects of a variety of systematic errors on estimates of 
the CMB polarization power spectra and have applied 
the method to a variety of possible errors. Let us begin 
by summarizing these results in a more compact form. 

To illustrate the relative magnitudes of the various 
sources of error, let us consider a fiducial set of exper- 
imental parameters. Let us assume that the true power 
spectra in the range of multipoles probed by our experi- 
ment are in the ratio 

C"^ : : = 300^ : 300 : 1, (7.1) 

roughly typical for subdegree-scale experiments. Fur- 
thermore, let us assume a fiducial value of 

s2=0.02, (7.2) 

which corresponds roughly to a baseline formed by a pair 
of antennas separated by three times the antenna diam- 
eter. 

Having chosen these fiducial values we can work out 
the effect of any particular error source. For instance, 
consider the effect of gain errors on a linear experiment 
aiming to measure B polarization. The leading contribu- 
tion to the error is the one that couples EE to B, with 

^B.EE,., = = 0-071. (7.3) 



The effect on the measurement of C is 

^=^.,..,.2 = 2172. (7.4) 

Say for instance that we wish systematic errors to have 
at most a 10% effect. Then 2I72 < 0.1 or 72 < 5 x 
10^'^. Of course 72 here represents the r.m.s. value of an 
unknown gain fluctuation, so this should be interpreted 
as an estimate of the level to which gain fluctuations must 
be understood. 

Tabled summarizes the results of such calculations for 
the various errors considered in this paper. A horizontal 
line separates instrument from sky errors. In each case, 
the dominant term listed is the one that involves the 
largest input power spectra. In cases where is small, 
an error term that is lower in the hierarchy may be of 
comparable significance to the dominant term. The table 
therefore lists a second contribution to each error where 
appropriate. This second contribution has k more weakly 
dependent on than the dominant contribution, so for 
large antenna separation it may be the more important 
term (although for the fiducial parameters adopted here 
it never is). In the cases of coupling errors and cross- 
polarization in an £'-mode measurement, the dominant 
term is independent of , so there is no need to consider 
a second term. 

In all entries in the table, the coefficients are averaged 
over a and calculated with the leading-order term in an 
expansion in s^, As figures indicate, the latter ap- 
proximation is excellent. 

In all cases, the error parameters should be taken as 
r.m.s. residuals after known errors have been removed. 
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For instance, as we noted in the previous section, sky 
curvature can induce cross-polarization characterized by 
/i = icT^. Presumably that effect would be known and 
accounted for; the parameter /i in Table represents an 
unknown and hence unmodeled additional component. 

Not surprisingly, the coupling parameters e and cross- 
polarization fi are of the greatest concern, since they cou- 
ple the temperature power spectrum to polarization mea- 
surements, in particular, if we want SC^ /C^ to be, say 
at most 10%, then these parameters must be e, /i < 10^^. 

Recall that for a linear experiment the coupling pa- 
rameters can be used to describe errors in the alignment 
of the polarizers, so a B mode experiment would require 
alignment with a precision ~ 10"'' radians or ~ 0.3'. For 
the E power spectrum, on the other hand, the required 
tolerance is about 0.3°. 

For pointing and beam shape errors, circular exper- 
iments have an advantage over linear experiments, be- 
cause errors that differ between measurement of Vq and 
Vu (parameterized by (5_,C-) are absent. Gain errors, 
on the other hand, are worse in a circular experiment. 

All of the errors in Table HI are expressed as couplings 
between band powers. In the case of instrument errors, 
we have seen that the visibility errors can be expressed 
as linear combinations of the visibilities themselves. In 
other words, the Fourier-space window functions associ- 
ated with the errors have exactly the same shape as the 
visibilities themselves. In the case of beam errors, this 
is not strictly true: 5Vq, for instance, has a different 
window function from Vq. However, for all of the er- 
rors considered in this paper, differences in Fourier space 
sensitivity introduced by the errors are relatively small: 
in all cases, the errors sample regions of Fourier space 
centered near k = 27ru with widths Au ^ cr^^, just as 
the visibilities themselves do. In short, the errors do not 
couple greatly different angular scales to each other. This 
contrasts with single-dish imaging experiments, in which 
scale-scale coupling induced by systematic errors is an 
important consideration |l7j . 

The results above were calculated using a simple and 
relatively conservative model for propagating errors from 
visibilities to power spectrum estimates. In a real data 
set, each resolution element in the Fourier plane would 
be sampled by multiple visibilities rather than just one 
pair. If we can assume that the errors in all of these vis- 
ibilities are independent of each other and have "nice" 
probability distributions (particularly that the errors are 
centered on zero), then the estimates should be reduced 
by a factor of -v/lv where N is the number of independent 
visibility pairs (Vq, Vu) per resolution element. However, 
since systematic errors often do not have nice statistical 
properties, a more conservative approach may be war- 
ranted. Even if errors do not need to be removed to the 
levels indicated here, it seems safe to say that their prop- 
erties need to be studied down to this level of precision 
in order to have confidence in the results. 

Although there is expected to be no cosmological cir- 
cular polarization in the CMB, it is worthwhile to con- 



sider the effects of circular polarization in the context of 
systematic errors. On the one hand, various errors can 
couple any intrinsic circular polarization that does exist 
(e.g., from foregrounds) into the linear polarization chan- 
nels, resulting in spurious E and B signals. On a more 
positive note, assuming that there is no intrinsic circular 
polarization, monitoring the circular polarization visibil- 
ities Vv may provide a way to assess systematic errors. 
In particular, in a linear experiment coupling errors (in- 
cluding polarizer misalignments) lead to a contribution 
to Vv that is correlated with the temperature anisotropy. 
Considering the level of control of these errors that is re- 
quired in a B mode experiment, such a diagnostic may 
prove quite useful. 
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APPENDIX 

In Sections IIV CI and IVII we made use of the following 
fact: Let v be a complex Gaussian random vector with 
mean zero and covariance matrix 

M=(v-v^), (A.l) 

and let A be an arbitrary hermitian matrix. Let q be the 
quadratic form 

g = vt • A • V. (A.2) 

Then the mean-square value of q is 

(q^) = Tr[(A • M)2] + [Tr(A • M)]^. (A.3) 

This Appendix, which no one will ever read, provides a 
proof of this fact. 

First, note that we can always reduce the problem to 
an equivalent one in which M is the identity matrix. To 
see this, let Q be a matrix such that M = Q • (e.g., 
by Cholesky decomposition). Let v' = Q^'v and A' 
Qt • A • Q. Then g = v't • A' • v' and (v' • v't) is the 
identity matrix. We will assume that this transformation 
has been made and drop the primes. 

Now diagonalize the hermitian matrix A: 

A = Rt • A • R, (A.4) 

where A is diagonal with real entries A, , and R is unitary. 
Let v' — R-v. The covariance matrix of v' is the identity 
matrix: 

(v' • v't) = R • (vvt) . Rt = R . Rt = 1. (A.5) 
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We have 



9 = v't.A.v' = ^A,K| 



and therefore 



(A.6) 



(A.7) 



Each of the quantities v'^ is an independent complex 
Gaussian random variable with mean zero and variance 



one, so 



[2 i{i=j 



(A.i 



(For real numbers, the i = j case would be 3 rather than 
2.) 

Writing this as 1 + dij , we conclude that 

= AiAj + XiXjSij (A.9a) 



[TV(A)]2 + TVCA^). 



(A.9c) 



Since traces arc unchanged under similarity transforma- 
tions, Tr(A) = Tr(A) and Tr(A2) = Tr(A2). We have 
thus established the desired result. 
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